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Abstract 

We have investigated a general structure of the ground-state wave 
function for the Schrodinger equation for identical interacting 
particles (bosons or fermions) confined in a harmonic anisotropic 
trap in the limit of large A^. It is shown that the ground-state 
wave function can be written in a separable form. As an exam- 
ple of its applications, this form is used to obtain the ground- 
state wave function describing collective dynamics for A^ trapped 
bosons interacting via contact forces. 
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The structure of the ground-state wave function for a many-body sys- 
tem is very important for theoretical understanding of recently observed 
Bose-Einstein condensation (BEC) [1] (the theoretical aspects of the BEC 
are discussed in recent reviews [2]) and other many body problems. The 
Ginzburg-Pitaevskii-Gross (GPG) equation [3] is most widely used to de- 
scribe the experimental results for the BEC. Recently, an alternative method 
of equivalent linear two-body (ELTB) equations for many body systems has 
been developed based on the variational principle [4,5]. In this paper, we 
consider identical particles (bosons or fermions) confined in a harmonic 
anisotropic trap. We show that in the case of large N the ground-state wave 
function can be written in separable form as 



^'(fi, f2, ...rW) = (pix, y, z) ■ x{^, <J), 
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Q is a set of (3N - 3) angular variables, and cr is a set of spin variables. 

We start from a generalization of the hyperspherical expansion of the 
wave function for the Hamiltonian 



^2 AT 1 ^ 

^''^ i=i ^ i=i i<j 

in the form [4,6] 



(3) 



v^(ri, ...r^) = Y: ^[k]{x, y, z)Y^k]{^^ , , Of, a), 

[K] 



(4) 



where Y^k]{^^ , , l^f , a) = Y^Xki^^ ^ ^ ' ^) ^^e combination of 
the hyperspherical harmonics, ¥^^(0;^), Y^^(f2^), and y^^(n^), with func- 
tions of spin variables a, which is symmetric or antisymmetric with respect to 
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permutations of particles for bosons or fermions respectively. [K] represents 
a set of numbers [K^, u^, Ky,i/y, Kz, Vz]- 

The hyperspherical harmonics Y^^ {Q,^ ) , Y^^{Vty ) , and Y^^ (Q^ ) are eigen- 
functions of the hyperspherical angular parts of the Laplace operators Z^^i 
Eili and Eili respectively. 

The Laplace operators are defined by 
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^ay/ ^dy dy y^ ^ 



and 

^ ^2 Id (9 1 

The hyperspherical angles 9f,92,...0%_i,9i,92,...0'%_i,9l,92,...0f^_i can 
be chosen in such a way that the hyperspherical angular parts of the Laplace 
operators Aqn satisfy the recursion relation [7] 

"^^^ = sin^-^]^_,a|^^''"''"^^-^a|:;^ + ih^V-- (6) 

with u — x,y, or z. 

Functions ^[K]ix,y, z) satisfy equations 



J2kmK']^[K']{x,y,z) = E^iK]{x,y,z), (7) 

[K'] 
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where 

h[K]lK'] = ^K^K'jKyK'^^K^K'^^v^v'Juyu'yKv'^ — i^ ^ ^ 5^2 ) 

2 2 2 2^^ ^ {N-l + 2K,){N-3 + 2K,) 

+ +U^yy +^^^) + ^( ^2 

{N -l + 2Ky){N -3 + 2Ky) {N - 1 + 2K,){N - 3 + 2K,) 
+V[K][K'](x,y,z), 

(8) 

with 

V[K][K'](x, y, z) =< K^, u^, Ky, Uy, K„ \ ^ l^i„t(r— r^) | K'^, v'^, K'y, u'y, K'^, v'^ > 

(9) 

We write ^[K]{x,y, z) in the form of a Laplace integral 

^[K]{x, y,z)^ j f[K]{oix, oiy, a^)(f)x{x, ax)(t)y{y, ay)(j)z{z, az)daxdaydaz, 

(10) 

where 



and u — {uJxUJyUzY^^- 

The Hill- Wheeler type equations [8,9] are obtained by requiring that en- 
ergy of the system is stationary with respect to the functions fiK]{(^x, otyi O-z) 
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^ / daxdaydaJ[K]{ax,(yy,az)[H[K][K']{oixOiyaz,a'^aya'^) 

(12) 

-S[K][K']S(axayaz, a'^a'ya'^)E\ = 0, 

where 

H[K][K']{otxOiyOiz,a'^oCya'^) =< (t)a:{x,aa)(t)y{y,ay)(t)z{z,az)Y]^K] 

(13) 

X I I (t)x{x,a'^)(t)y{y,a'y)(l)z{z,a'^)Y[K'\ >, 

and 

^(aajaya^,^^^^^^) =< (px{x,ax)4>y{y,ay)4>z{z,az) \ 4>x{x,a'^)4>y{y,a'y)4>z{z,a^) > 

(14) 

In order to solve the Hill- Wheeler type equations (12), we assume that 
the integral in Eq. (10) can be replaced by sum 

oo 

<^>^K]{x,y,z) = c^S^x{x,al)(py{y,al)(pz{z,a';), (15) 

i,j,k=l 

where c|^' are solutions of the following equations 

Yl [H[K][K']{.Oilala'l,aia^ya'i) - 5[K][;^']'S(Q;>^a^ ^'a^ )-E]cl5l' = 0. 

i',3',k' 

[K'] 

(16) 

For the case of large N, the overlap, Eq. (14), 

K . y . y L((«^)2 + (4')2)((«.)2 + («.)2)((«.)2 + («.')2)J 

(17) 
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reduces to the Kronecker deltas 

S{aiaia';, o^a^a'; ) = Sa^SjfSkk' (18) 

Since the ratio 

is a much more slowly varying function of a compared to ^'(a^a^a^, a^^a^'ya'^') 
in almost all cases [10], we have for the case of large N 

H[K][K'] {ala^a';, a^'a^'af ) = H^khk'] {(ix, dy, (iz)Sii'Sjf6kk', (19) 

(see Appendix for the case of N identical particles interacting via contact 
repulsive force). 

Substitution of Eq. (19) into Eq. (16) gives 

^[K]{x,y,z) = C[K](px{x,d:c)<Py{y^<^y)M^:(^z): (20) 

where C[k] are solutions of the following equations 

^[H[K][K']idx, dy, dz) - S[K][K']E]C[K'], (21) 
[K'] 

and parameters dx,dy, and are solutions of 

dE _ dE _ dE _^ 
ddx ddy ddz 

Substitution of Eq. (20) into Eq. (4) yields ^(ri,r2, ...rjv) given by Eq. (1) 
with 

<t>{x, y, z) = 0^(x, dx)(t)y{y, dy)(f)z{z, d^), 



and 



[K] 



We now consider identical particles confined in an anisotropic harmonic 
trap and interacting via contact force 

VintKTi - Tj) = b[ri - Tj), (23) 

with positive scattering length a > 0. Using factorization (1) we have 

^^iN,:m^^2^ + ^ + ^) + ^]<Pi^, y, - m^, y, z), 

(24) 

where =< x | Aqat \ x> / < x \ X> with t = (x, y, z) and 
_ a^^7V(7V-l)^ r(iV/2) 



v/2^m T((Ar-l)/2)' 

In the large A?" limit, parameters Cx,Cy,Cz, and c are expected to be slowly 
varying functions of N. For N identical bosonic atoms with large N, an 
essentially exact expression for the ground state energy can be obtained by 
neglecting the kinetic energy term in the GPG equation [2,3] (this is called 
"Thomas-Fermi approximation" [11]). From comparison of the ground state 
solution of Eq. (24) with the Thomas-Fcrmi approximation [11], we can fix 
unknown parameters and find the ground-state solution of Eq. (24) as 

H^, y, z) = i)x{x)il)y{y)il)^{z), (25) 

E = — - — n^/^ (26) 
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with 

ij,{x) = Ax^^-'^/^cxp[-mi:;{x/af/{2h)], 

i:y{y) = exp[-mu{y/pr/{2h)], (27) 

ilj,{z) = Az^''-^^/^exp[-mu{z/^f/{2h)], 

where A = ^2/T{N/2){mCj /{a^h))^/^, a = n}'^Lb/uj^, (3 = n}'^Cj/ujy, 7 = 
n^f^uj /ujz-i = {'j^x'^y'^zY^^ ■I'l^ — nCjU = 2^Jujm/ {2nh)Na and 

c=(^r/^^v^«^0.82. (28) 

Eqs. (25-28) give the exact ground-state solution of Eq. (24) for large N . 
Thus we have found an analytical solution for the ground-state wave function 
describing collective dynamics in variables (x,y,z) in the large N limit. 

We note that the slope of the Thomas- Fermi wave function becomes in- 
finity at the surface, leading to logarithmic singularity in the kinetic energy. 
Hence it is necessary to modify the Thomas-Fermi wave function near the 
surface [12-14]. In contrast, we do not have such problems for our solution, 
Eq. (25-28). 

It is also interesting to compare our results with the ELTB method [4,5]. 
For this situation (contact force, Eq. (23) and large N limit), the ELTB 
method corresponds to c^/^ = 1. It shows that the ELTB method is a very 
good approximation with relative error of about 8% for parameter c^l^ . 

In summary, we have investigated the general structure of the ground- 
state solution of the Schrodinger equation for N identical interacting particles 
(bosons or fermions) confined in a harmonic anisotropic trap in the large N 
limit. The main results and conclusions are as follows 

(i) It has been shown that in the case of large N the ground-state wave 
function can be written in separable form, Eq. (1). 
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(ii) Using this form, we have found an analytical solution for the ground- 
state wave function, Eqs. (25-29), describing collective dynamics in collective 
variables (x,y,z) for N trapped bosons interacting via contact repulsive forces. 

(iii) Our results can be used for checking various approximations (both 
existing and future) made for the Schrodinger equation describing N identical 
interacting particles (bosons or fermions) confined in a harmonic anisotropic 
trap. 
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Appendix 



To prove Eq. (19) we consider the contact potential case 

Vi„t{ri - rj, a) = S{ri - rj)r){a), {A.l) 

where rj depends on spin variables. Using Eq. (A.l) we can rewrite Eq. (9) 
as 

N(N - 1) 

VlK]lK']{x,y,z) = -flK][K'] — , (A2) 

where 7[i<:][i<r'] does not depend on x, y, z. 

Substitution of Eq. (A.2) into Eq. (13) gives 



T{N/2) 



(^.3) 

where Pt — <^t/<^ for t — x,y, or z. 

For large N, S{al.a^ya'^, al^a^yaf), Eq. (14), reduces to the Kronecker 
deltas 6ii'6jj'6kk', and hence from Eq. (A. 3) we obtain Eq. (19). 
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